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1. Introduction 

The term "microhalo" loosely refers to a generic dark matter, gravitationally-bound 
object, whose mass is typically less than that of the Sun. Objects of sub-solar mass 
scales are predicted to be the first objects formed in the high-redshift Universe in 
theories where the dark matter particle is "cold" (such as ACDM). In this short article 
I summarize the present understanding of microhalos. In Sec. [2] I review the theoretical 
motivation behind the presence of microhalos, in Sec.[3]I present the results of numerical 
simulations that studied their formation. Sec. H] is an overview of the present state of 
understanding the survival of these objects, while in Sec. [5] I discuss the procedure of 
how to characterize their physical properties. I present their connection to direct and 
indirect detection experiments in Sees. [6] & [3, and I conclude in Sec. EJ Throughout 
the paper, the assumed cosmology is flat ACDM, with VL m h 2 = 0.1358, VL^h 2 = 0.02267, 
h = 0.705, cr 8 = 0.812, and a spectral index of n s = 0.96 pp. 

2. Kinetic decoupling and microhalo scales 

In the standard cold dark matter cosmological model, the energy density in the Universe 
is balanced by approximately 4% baryonic matter, 23% cold dark matter and 73% 
dark energy. Cold dark matter refers to gravitationally interacting matter that is 
hypothesized to be in the form of a yet-to-be-discovered particle. From the particle 
physics aspect, new particles arise in essentially all extensions to the standard model of 
particle physics, and usually the lowest mass particle is stable due to a new symmetry. 
Dark matter particle candidates include (among many others), neutralinos, the lightest 
Kaluza-Klein particle, axions, and sterile neutrinos. Even though there is no a priori 
physical reason as to one candidate being more favorable than another, neutralinos 
are attractive and well-studied because they are experimentally accessible at present. 
Neutralinos arise in supersymmetric extensions to the standard model of particle physics 
[21 [3], and they are part of a generic class of dark matter candidates, called WIMPs 
(Weakly Interacting Massive Particles). 

In the early Universe, conditions are such that a WIMP is in chemical, thermal 
and kinetic equilibrium. As the rates of these interactions are diluted due to the 
expansion of the Universe, the WIMP falls out of equilibrium and decouples. The 
decoupling temperature is determined by the scattering cross section that is responsible 
for the equilibrium condition the particle decouples from, while the velocity distribution 
function of the particle is set by its mass. The temperature at which a WIMP is 
kinetically decoupled, i.e., when momentum-changing interactions cease to be effective, 
is called the kinetic decoupling temperature. For supersymmetric dark matter, such as 
the neutralino, the kinetic decoupling temperature is T d ~ [10 — 1000] MeV [HG3E]. The 
free-streaming of particles after kinematic decoupling tends to smooth out fluctuations 
on small scales [H HI [8j [9] . This physical effect leads to a cutoff in the dark matter power 
spectrum. 
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Figure 1. Left: An analytic example which shows the damping factor of the power 
spectrum [5J[H]- 

The cutoff of the power spectrum includes an imprint of the acoustic oscillations of 
the cosmic radiation fluid [10] . In general the transfer function of the dark matter 
density perturbation amplitude is obtained by solving the Boltzman equation, but 
under certain approximations, it is possible to be derived analytically [H [9]. Fig. CD 
shows an example of the analytic calculation of the cutoff of the power spectrum, 
derived from [8j [9]. What is shown is the damping factor D(k) that arises for a 
super symmetric WIMP dark matter particle with a decoupling scale of fed ~ 40 pc -1 , 
and a free-streaming scale of k ls ~ 0.94 pc -1 (corresponding to a kinetic decoupling 
temperature of Td ~ 20MeV). In general, the cutoff scale is related to the kinetic 
decoupling temperature as M cut « lO _4 M (T d /lO MeV) -3 [10]. Recently, [11] presented 
the posterior of the cutoff scale from a Bayesian analysis of the presently-allowed 
supersymmetric parameter space. 

3. Microhalos in numerical simulations 

3.1. Formation of microhalos 

The first attempt to numerically simulate the formation of microhalos was done in 2005 
by the group at the University of Zurich [12]. In that work, a small region of the Universe 
was simulated using a multi-scale technique [13] in order to achieve the required spatial 
resolution needed to resolve objects below the cutoff-scale of the power spectrum. The 
input transfer function was taken from [8]; corresponding to a supersymmetric dark 
matter candidate of mass ~ lOOGeV. This particular choice leads to an exponential 
cut-off at ~ 1O _6 M . The simulation was evolved from a redshift of z — 350 to a 
redshift of z = 26, when the simulated region itself approached non-linearity. In the 
small, high-resolution volume of roughly average density, the first non-linear structures 
with mass near the cutoff scale of ~ 10 _6 M Q are observed to form by redshift z ~ 60. 




The detection of sub-solar-mass dark matter halos 



4 



Their physical size is few x pc, and they are described by a single power-law density 
profile. Even though halos of mass in the range [10~ 6 — 10°]M Q are formed roughly at the 
same time, they all share the same concentrations at z = 26. As expected, these objects 
contained no substructure, and no objects with smaller mass were formed, confirming 
the input shape of the power spectrum. Even though the statistics in this first simulation 
were low, the deduced mass function suggested a power law and normalization consistent 
with the halo mass function derived from numerous orders of magnitude larger scales, 
namely, dN/dlnM ~ M _1 . This finding is an important piece of information when it 
comes to studies aimed at the long-term survival of microhalos in the potential well of 
a present-day galaxy, such as the Milky Way. 

In a subsequent paper [14J, the numerical simulation of the formation and evolution 
of microhalos was studied by using 64 million dark matter particles within a comoving 
volume box of 3 kpc on the side (implying a particle mass of 9.8 x 1O _1O M ). The initial 
redshift was z = 456, and the initial conditions were similar to the simulation discussed 
in [12] (see previous paragraph). By the end redshift of the simulation, z ~ 75, the 
volume contained almost 2000 virialized dark matter sub-solar mass halos. The most 
massive halo is ~ 1O~ 2 M , and corresponds to a 3.5a fluctuation in the density field. For 
a comparison of the evolution and substructure content of this simulation, the authors 
performed a separate simulation, of similar dynamic range, of a cluster-size halo at 
z = 0. Comparison between the substructure mass functions in the two simulations 
showed confirmation of the earlier result found in [12], namely a power law behavior 
that is self-similar over more than 10 orders of magnitude in subhalo mass. 

It should be emphasized that the simulation presented in [12] was performed in 
a region of average density, thus probing the formation of "field" sub-solar mass dark 
matter halos. On the other hand, the high-resolution simulation presented in [13] was 
centered around the highest density peak of the initial 3 kpc volume. This allowed the 
study of the assembly and structure of an object with mass higher than the cutoff scale 
of the power spectrum, while resolving scales down to the cutoff scale. In addition, 
it provided the first estimate of the survival rate of dark matter microhalos, a very 
important result which relies on the non-linear evolution of structure formation. 

3.2. Challenges in numerical simulations 

The numerical simulation of the formation and survival of microhalos is a very difficult 
task. The difficulties arise from two main sources. First, by the fact that the power 
spectrum P(k) of sub-solar mass scales asymptotically approaches the k~ 3 behavior, 
thus leading to the formation of multiple scales roughly at similar epochs, and second, 
by uncertainties in the actual implementation of numerical techniques in numerical 
simulations [15] . 

The cutoff in the power spectrum that enters the initial conditions in the numerical 
simulations does not fully suppress the formation of structure at scales smaller than the 
cutoff scale. Even though the fraction of substructure in objects below the cutoff scale 
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is small, their mere presence deserves further study in order to determine whether this 
effect is an artifact of the simulation techniques employed, or whether it is a physical 
effect. It should be noted that such effects have also been found in simulations of Warm 
Dark Matter (WDM), where the dark matter particle possesses a non-zero thermal 
velocity [EJ [T71 [TBI [T9 J . A possible numerical source that leads to this artificial effect 
is the implementation of grid initial conditions [20J. If not implemented at an early 
enough starting redshift, the initial Zel'dovich displacement can be a large fraction of the 
inter-particle separation [21]. This can lead to artificial fragmentation along filaments, 
however the effect has also been observed when the use of glass initial conditions [22] is 
implemented. The simulations of [H] also show that the formation of sub-cutoff scale 
structure is not limited to a preferred separation along a grid axis, suggesting that it 
may perhaps be a real effect. 

Another source that can introduce bias in the measurement of a mass function is 
the starting redshift of simulations. As shown in [21j, the requirement that the Fourier 
modes within the simulation box are all linear is very naive. The mass function cannot 
be robustly determined at redshifts which are close to the redshift of first crossing (the 
redshift where halos form). It has also been pointed out that even perhaps the choice 
of a halo finder in the simulation can also lead to an uncertainty in the measurement of 
the mass function [HI [15], due to the fact that the density contrast in sub-solar mass 
halos is low due to similar formation times. Therefore, caution must be taken in the 
interpretation of these results, and a thorough convergence study on the mass function 
of sub-solar mass halos is of paramount importance, and in great need. 

4. Survival of microhalos in the Milky Way halo 

The survival (and thus existence) of microhalos in the present-day Milky Way halo has 
been a subject of debate. The approach typically taken in these studies has progressed in 
two fronts: simple analytical estimates, and through numerical simulations. In principle 
the two should agree, however given the vast orders of magnitude in dynamic range and 
highly complex non-linear effects, it has proven to be a difficult task. 

The question of surviving microhalos was first addressed analytically in [23]. A 
detailed analytical study of the effects of tidal destruction on microhalos due to rapid, 
early-time hierarchical merger phase, showed that only [0.1 — 0.5]% of formed microhalos 
would survive to the present epoch [23J. The initial claim that a large number of 
microhalos could be present in the solar neighborhood [12J was challenged only days 
later [21]. The argument presented was that strong impulses by individual stars in 
the disc can tidally disrupt microhalos, leading to a large number of tidal streams. 
Within days, a response [25] showed that the application of the impulse approximation 
requires microhalos to experience multiple stellar interactions and as a result the survival 
timescale is many Hubble times. 

One of the key advantages of numerical simulations in the study of the formation 
and abundance of microhalos is the ability to trace the history of each bound structure. 
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In doing so, the work of [TJ] found that substructure formed at early times is highly 
prone to tidal disruption due to the high central density of the host halo. In addition, 
most of the remaining microhalos that survived to the end-redshift of the simulation 
run are significantly tidally stripped. As expected, the amount of mass-loss experienced 
due to tidal effects depends on the history of the subhalo and in particular whether 
it experienced passages near the center of the host halo. However, most importantly, 
the study of [H] demonstrated that even though the density contrast of sub-solar mass 
scales is small, microhalos of mass 10 _6 M Q are able to survive in a host of mass at 
least 10 4 times more massive, thus addressing numerically the question of survival in 
the initial merging process. 

A separate study [26] looked at the energy input to microhalos due to stellar 
interactions by using analytic calculations, and then tested their results in a numerical 
simulation [27J EE]- The main result is that at large impact parameters, the energy 
input to microhalos is independent of mass. They show that the survival timescale 
of microhalos in the Milky Way is of order the age of the Milky Way, but decreases 
significantly as the microhalo mass increases. In an independent test of the impulse 
approximation |29j, it was shown that multiple stellar encounters can remove most of 
the mass from microhalos, however the very inner cores can survive to the present time. 
Nevertheless, the normalization and shape of the density distribution of dark matter 
in microhalos will be decreased, with significant implications to direct and indirect 
detection experiments (see following sections). It is important to realize that it is 
currently not feasible to fully address the issue of survival, mass and spatial distribution 
of microhalos, due to the enormous dynamic range required and the limited computing 
power that technologically exists. Studies such as the ones mentioned here are invaluable 
in shedding light on what might happen in ideal situations, and therefore should be 
further explored. 

5. Characterizing the physical properties of microhalos 

In order to assess the detectability of microhalos in any scheme (e.g. direct or indirect 
detection experiments), it is first necessary to determine the distribution of dark matter 
within the virialized region that comprises the microhalo. Typically, this can be done 
by assuming a profile, and a normalization of that profile. 

The numerical work of [T2], [TJ] finds that microhalos at the redshift at which they 
were studied in the simulation are described by an NFW profile [30]. The form of this 
profile is p(r) = p s /f (1+f) 2 , where p s is a density normalization, called the characteristic 
density, r s is the scale radius, and f = r/r 8 . The two parameters needed to specify the 
profile are the characteristic density, and the scale radius. The numerical simulation 
findings suggest that r s is approximately equal to the virial radius i? vir of the microhalo, 
suggesting an NFW concentration parameter (c v = i? v ; r /r s ) at that redshift which is of 
order one (more specifically c v ~ [2 — 4.5] [31]). For field microhalos, the concentration 
parameter grows as (1 + z) -1 , indicating that "field" subhalos have concentrations in 
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Figure 2. Left: The rms fluctuation on sub-solar mass scales for the WMAP5 
cosmology. Right: The concentration parameter of "field" sub-solar mass halos at 
2 = based on the model of 33J . This is obtained from an extrapolation of the c v (M) 
relationship found on galactic scales, all the way down to microhalo scales. It does 
not represent the c v (M) relationship of Milky Way progenitors that were of microhalo 
scale at the time of their formation. 

the range of [50 — 100] for a formation redshift around z ~ 20 (see below). Obviously, 
the profile of a present-day microhalo in the Milky Way has most likely evolved from 
this initial distribution of dark matter, to a distribution which is the result of tidal 
interactions throughout the course of structure formation. As I alluded to in Sees [3] & 
HI these processes are very difficult to accurately model. Nevertheless, it is a reasonable 
assumption to consider the initial profile at the redshift of formation as a representation 
of a maximal distribution, and keep in mind that any experimental ramifications are 
going to actually be upper limits and not absolute values. 

In order to normalize the dark matter density profile of microhalos we need to first 
determine their formation time, which can be obtained by knowledge of the CDM power 
spectrum, and the cosmological growth factor. 

The normalization of a the profile of a dark matter halo can be obtained by assuming 
that mean density of the halo is A vir times the mean matter density of Universe pyi, 
with A vir being the virial overdensity (see e.g. [32]). The mass of the halo is then 
M = 47ipM^vh(z)Rl ir /3. By using this definition of the mass of a dark matter halo, 
the normalization of the NFW profile is then simply obtained via p s = puA v i r cl/3f(c v ), 
where f(x) = ln(l + x) — x/(l + x). The rarity of a particular mass scale at a particular 
redshift can be obtained from the power spectrum of fluctuations. In general, if we 
define the power in each logarithmic interval in k as A 2 (k) ~ k 3 P(k), where P(k) is 
the CDM power spectrum, then the mean square fluctuations on a scale M is given by 
a 2 (M) = J A 2 (k)\W(k; M)\ 2 d\n k. Here, W(k; M) is the Fourier transform of a real- 
space spherical top-hat window that contains mass M. The left panel in Figure [2] shows 
the root mean square fluctuations on microhalo scales, derived by assuming the WMAP5 
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Figure 3. Left : The formation redshift (according to linear theory) as a function of 
halo mass. The curves correspond to 1 (solid red), 2 (short-dash blue) and 3 (long-dash 
green) a fluctuations in the density field. Right : The NFW characteristic density of 
halos collapsing out of 1, 2, and 3 a peaks. In all cases, the underlying cosmological 
model is taken from the combined analysis of WMAP5 with distance measurements 
from the Type la supernovae, and the Baryon Acoustic Oscillations [T]. 



cosmological parameters pQ. Note that for illustration purposes, the primordial power 
spectrum in this example does not contain the cutoff shown in Fig. [TJ Knowledge of the 
underlying cosmological model (i.e., the values of the cosmological parameters), allows 
us to compute the growth function of fluctuations D(z). Then, the redshift of formation 
of a mass scale M is defined as the redshift z where a(M)D(z) = 5 C , where 5 C = 1.686 
is the value of the characteristic density for collapse (derived from linear theory). The 
rarity of a particular peak is typically defined to be in units of the standard deviation 
of the smoothed density distribution, so that v = 5 c /a(M)D(z). 

Note that the rarity of a particular fluctuation in the density field is directly related 
to the abundance of that scale in the present Milky Way halo - a key quantity in any 
dark matter experimental effort. For example, a smaller fraction of the total Milky Way 
mass originates from 3a peaks in the density field as those are rarer than e.g. la peaks. 
However, the 3a peaks are formed earlier, thus they will be denser. As such, they are 
not as susceptiple to tidal disruption as the la peaks. The abundance of microhalos 
in the present epoch will be determined by a balance of these competing effects and a 
full description of the microhalo population of the Milky Way must include all of these 
factors. A first step in that direction was the study of [34], where it was found that 
high-z/ progenitors are found predominantly found near the centers of dark matter halos 
today. For example, according to |34j . the median progenitor at the solar neighborhood 
of the Milky Way corresponds to a 2a peak. 

An additional piece of information which is of importance when it comes to dark 
matter detection experiments, is the tidal radii of survived microhalos in the Milky Way. 
A reasonable approach would be to assume that the truncation radius at the present 
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epoch can be estimated by finding the radius of the microhalo where the density is equal 
to the density at the solar radius (35] . This assumption results in the truncation radius 
being equal to the virial radius for the microhalos halos considered here, reflecting the 
dense state of the Universe at the time they were formed. Clearly, the assumption here 
is that the halo profile normalization is intact from the time it was formed until the 
present epoch. Again, this points out at the lack of knowledge of the mass function and 
history (thus properties) of microhalos in the solar neighborhood. But given the high 
initial densities reflected in microhalos, and the steep power law of the density profile, 
it seems likely that the tidal radius will not be very different from their virial radius set 
at the time of collapse. 

A more general note regarding the physical description of the dark matter profiles 
of microhalos is with regards to the concentration parameter. It is dangerous to assign 
concentration parameters that would correspond to present-epoch "field halos" (i.e., a 
microhalo formed at high redshift and evolved uninterrupted to the present day). Field 
halos grow hierarchically by the accretion of smaller halos, and by smooth accretion of 
dark matter. It is not clear exactly how this process occurs at very high redshifts on 
microhalo scales, as the flatness of the power spectrum implies the simultaneous collapse 
of many different scales above the cutoff in the power spectrum. The maximum value of 
the concentration parameter for a field microhalo can be obtained by extrapolating the 
concentration-mass relationship deduced from numerical simulations [33, 36J (see Fig. [2]). 
Even though this is a daring extrapolation over numerous orders of magnitude, it can be 
considered as an upper limit for field microhalos, as the concentration-mass relationship 
stems from the shape of the power spectrum. The power spectrum flattens out at 
large k, and so is c v with respect to mass [33] (see Fig. [2]). The danger from assuming 
concentration parameters of field halos in microhalo studies near the solar neighborhood 
arises due to the fact that the density profile normalization is rather sensitive to the value 
of the concentration parameter p s ~ c^, therefore assuming concentration parameters of 
order 100, can result in a very large overestimate of the normalization of the profile. 

To summarize, given the uncertainty in the evolution of the profile of microhalos, 
the only reasonable approach is to determine their physical properties at the time they 
were formed, and consider those as upper limits to what they could be at the present 
epoch. Figure [3] &H] shows the formation redshift, characteristic density and virial radius 
of sub-solar mass dark matter halos. It is of vital importance to emphasize that any 
treatment of annihilation signals from microhalos (see e.g. 7-rays in the next section) 
must be consistent with a description of the physical properties of microhalos, which are 
set by the underlying structure formation paradigm. Arbitrary choices of concentration 
parameters (which affect the normalization of the profile), or inconsistent choices of 
formation redshift will only result in an estimate for a particular toy microhalo model. 
Such an approach cannot assess the structure formation relevance, and/or distribution 
function of expected signals and must treated with caution. 
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Figure 4. Left: The virial radius of microhalos at the redshift of collapse. Right: The 
annihilation luminosity of a sub-solar mass dark matter halo as a function of its mass. 
If the halo has lost 99% of its mass due to tidal interactions with other halos and/or 
baryons (stars, disk, clouds, etc.), then the luminosity should be decreased by a factor 
of 30% from the values shown. Line types and colors as in FigO 



6. Detection of microhalos using dark matter annihilation products 

It might be possible to detect the presence of microhalos in the Milky Way halo 
by searching for the annihilation products of the WIMP dark matter particle, and 
specifically for 7-rays. Under the supersymmetric WIMP scenario, the dark matter 
particle couples to quarks and leptons at the tree level, and to photons and the Z° 
gauge boson via one-loop diagrams. The search of 7-rays from the annihilation of the 
dark matter particle is attractive for two reasons: 1) a swath of 7-ray experiments 
are operational (e.g. Fermi/GLAST [37J, VERITAS [38J, H.E.S.S. [39], MAGIC 00], 
CANGAROO [H]), and 2) the shape of the emitted spectrum can be estimated (given 
a supersymmetric dark matter candidate). 

In general, the luminosity of a subhalo as a function of energy is given by 
L(E) = S JP(E)dE } where S = J p 2 (r)d 3 r (see Fig.HD, and P(E) = (dN '/ dE) (av) / 'M 2 . 
Here, dN/dE is the energy spectrum of the annihilation products under study (e.g. 
photons, positrons, neutrinos, etc.), (av) is the annihilation cross section, and M x is the 
mass of the dark matter particle. Maximal values of / P(E)dE in supersymmetric 
theories with a valid dark matter candidate (e.g. the neutralino) yield a value of 
< 10 -28 cm 3 GeV -2 for a threshold energy of ~ 1 GeV. 

6.1. Microhalos as 'y-ray sources 

In this subjection I discuss the possible detection of microhalos with 7-rays, though 
the formalism can be well-extended to other annihilation products, such as positrons 
and neutrinos. In order for a sub-solar mass halo at a distance d and flux on Earth 
$ = L/And 2 to be detected above a detector threshold $0, it must be located at a 
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distance d < \J$/Qq. Such an example is shown in Fig. [5j for a supersymmetric 
WIMP dark matter candidate with P < 10~ 28 cm 3 GeV -2 , and a detector threshold 
of $o = 10~ 9 cm -2 s _1 . It should be noted that for subhalo mass functions that scale 
as dN/d\nMdV ~ M -1 , the typical visibility distance between microhalos of mass M 
scales as (r) ~ M 1 / 3 . However, as the luminosity of a halo goes as L ~ p 2 r 3 ~ M, 
the "visibility" distance of microhalos is proportional to <i max ~ M~ 1//2 , i.e., it is 
"easier" to detect more massive halos at larger distance [121 H3], than small halos 
nearby. This scalings will be altered if the dark matter particle's annihilation cross 
section is Sommerfeld enhanced [lU US, HEJ HE] . The increased annihilation rate due 
to low microhalo internal velocities leads to the possibility that small nearby microhalos 
will be more luminous than the canonical case where the annihilation cross section is 
independent of velocity. 

One approach to investigate the observability of a population of microhalos in the 
Milky Way is to Monte Carlo the distribution and properties of microhalos as followed 
in [19]. The key unknown in this approach (or any approach for that matter) is the 
abundance of microhalos as a function of their mass, as well as their central densities (as 
there is a a distribution of mass scales that collapse in each redshift interval - see Sec. |5]). 
The results of jl9] show that for a fixed concentration microhalos (i.e., isolated density 
peaks), it is highly unlikely that there should be any high signal-to-nose single-source 
detection of microhalos with the all-sky Fermi Gamma-ray Space Telescope (FGST, 
formerly GLAST). 

Another approach was followed in [50] . In that study, the uncertainty that stems 
from the lack of knowledge of the survival rate of microhalos was parametrized as a 
fixed contribution to the local dark matter density, so that a particular fluctuation peak 
(i.e., collapsed microhalo at a particular redshift) results in a certain abundance. Again, 
the main caveat here is the unknown survival rate. The result of [50] suggested that a 
large number of microhalos could potentially be visible, albeit the normalization of the 
number was directly deduced from the unknown density parametrization. 

One unique feature of the potential presence of survived microhalos in the Milky 
Way is the probability that they can be present in very short distances (sub-parsec) 
from the Solar system. Thus, it may be possible to observe a proper motion of the 7-ray 
signal [25l l50~t |5T] . The two crucial conditions for the observation of the proper motion 
of a subhalo is that it is bright enough to be detected at some high signal-to-noise level, 
and it must be close enough so that the proper motion exhibited over the life-time of an 
experiment is above the threshold angular resolution of the detector. The details of this 
potential observation were studied in [50] under the assumption that a certain fraction 
of the Milky Way mass at the solar radius is in microhalos which were all formed in a 
particular epoch. 

Any potential detection of proper motion implies the presence of a large number 
of "unresolved" microhalos at distances far beyond the flux-limited distance set by 
the sensitivity of a particular detector. The presence of a large number of microhalos 
below individual detection threshold would lead to a background radiation that would 
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Figure 5. The visibility distance where the emitted flux is greater than the threshold 
flux of Fermi/GLAST as a function of microhalo mass for different peaks of the density 
field (see text for details). Line types and colors as in Fig.[2l 

correlate with the angular and radial distribution of microhalos in the Milky Way. 
This was the subject of study in [52], where it was shown that the measured value 
of the 7-ray background measured by EGRET [53] is already placing stringent limits on 
the probability of proper motion detection (due to the implied low number density of 
microhalos in the Milky Way). Addressing the same issue, [51] showed that regardless of 
the assumptions of the inner density profile of microhalos, the EGRET measurement of 
the 7-ray flux from the Galactic center [55] (see [56], but also [ST] ) is strongly suggesting 
that the probability of detecting microhalos with measurable proper motion is negligible. 
Nevertheless, given the simplicity of this measurement, a search for the proper motion 
of 7-ray-only sources should be performed due to the valuable information contained in 
such a potential detection. 

7. Microhalos and the local Milky Way dark matter distribution 

The presence of substructure in dark matter halos leads to interesting consequences 
when it comes to direct detection experiments. In general, the rate of events in a direct 
detection experiment is proportional to T ~ p Q , where p & is the dark matter density 
in the solar neighborhood. Typically, the value of p is obtained from dynamical 
measurements of the structure of the Milky Way (e.g. [35]), with a canonical value 
of p & = 0.4GeVcm~ 3 . This value is obtained from averaging regions which are of 
order ~kpc, much larger than the regions probed in the duration of a direct detection 
experiment (sub-parsec scales), and it is susceptible to uncertainties that are introduced 
due to the unknown shape of the Milky Way halo (see e.g. [581 [59]). If there 
are fluctuations on sub-parsec scales (as what one might expect from the presence 
of microhalos), then direct detection experiments are prone to structure formation 
uncertainties. In this section I give an overview of the experimental implications of 
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Figure 6. The local dark matter density probability distribution function [60J. The 
peak of the distribution defines the dark matter density of the smooth distribution 
of dark matter, while the power low behaviour at higher densities arises from the 
presence of dark matter subhalos. This function is derived under the assumption that 
the substructure mass fraction in the Milky Way of all substructure in the mass range 
[l(r 10 Af Q - M MW ] is ~ 20%. 

fluctuations in the density field at the solar radius. 
7.1. Microhalos and the mean dark matter density 

It is expected that if a certain amount of mass in a dark matter halo is distributed in 
high-density regions, then the density of the smoothly-distributed dark matter will be 
lower than the case where all matter is distributed smoothly. The presence of sub-solar 
mass microhalos in the Milky Way raises one important question: what is the value of 
the smooth component if the substructure mass function extends down to microhalo 
scales, and what is the probability of the solar system being in an overdense region 
at any given epoch. The answer to both of these questions is directly related to the 
local flux enhancement in the annihilation signal (the so-called "boost factor"'), as the 
probability of being in a density enhanced region is inversely proportional to the boost 
factor. 

Addressing this question was the study presented in [60J. Under the assumption of 
hierarchical structure formation, it was shown that the derived local density probability 
distribution function (PDF) is positively skewed due to the presence of substructure, and 
that the peak of the distribution is always less than the canonical mean value of the local 
dark matter density in the solar neighborhood (see Fig. [6]). This conclusion is derived by 
a simple analytic model based on the scale-invariant nature of the hierarchical structure 
formation, as well as by investigating the implications of the distribution and profiles 
of dark matter halos as seen in numerical simulations. Both approaches showed similar 
results. These results were recently confirmed in numerical simulations by analyzing 
a suite of halos taken from the Aquarius Project [EU [62] and independendly, by an 
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analysis of the Via Lactea II Milky Way simulation [5H1 E3] • 

An interesting outcome of the approach of [60] is a derivation of the the local 
annihilation rate due to the granularity of the local dark matter halo. The local boost 
factor should not be confused with the global boost factor that is commonly referred 
in the litereature (see e.g. [Ml EI])- The least granular cases studied in [60J yield a 
boost factor which is roughly consistent with estimates based on the survival rate of 
microhalos [23J as well as direct simulation results, as shown in [63]. The local boost 
factor is ~ [1 — 5], with a weak dependence on the cutoff scale of the subhalo mass 
function. 

The introduced variance in the local dark matter density has significant implications 
to the interpretation of combined results from indirect and direct detection experiments, 
as the expected signal is proportional to different powers of the local particle density. 
The volume probed during a 3-year direct detection experiment is very small ( ~ 
few x 10~ 4 pc), while dynamical estimates of the local dark matter halo provide averages 
over much larger scales. The uncertainties implied by fluctuations in the local dark 
matter density due to the presence of substructure are manifested as uncertainties in the 
predicted rates in a dark matter detector |60j . If substructure in the form of microhalos 
is abundant in the solar neighborhood, then the local boost factor will be large 
(favoring indirect detection searches), however the probability of being in an overdense 
region will be low (disfavoring direct detection experiments). Fortuitously, present 
simulation results suggest that these effects (and similarly perhaps the uncertainties 
due to structure in phase space) are small, and a smooth dark matter halo is a safe 
assumption [611 [621 [591 E3], though the consequences for indirect detection searches 
may be more pronounced (e.g., the explanation of the postitron excess in the PAMELA 
[65] and ATIC [66] data by the presence of a nearby subhalo [67J). 

An interesting question would be to ask whether it is possible that a long-duration 
direct detection experiment will be able to map the local distribution of dark matter. 
If we assume the velocity of the Sun to be 220 km/s, then the distance sweept over 
~ 20 years, is x ~ few x 10~ 3 pc. An extrapolation of the subhalo mass function to 
microhalo masses (normalized so that 10% of the Milky Way halo is in objects with 
mass greater than 10 _5 Mmw) and assuming an NFW profile for the radial distribution 
of subhalos, results in a mean distance between microhalos which depends on mass as 
(r) pa 52M a3 pc. For microhalos of mass 10 _6 M Q , this is (r) pa 0.8 pc. As i C (r), 
it seems unlikely that a direct detection experiment will measure the transition from 
a smooth component to an overdensity of dark matter. Coupled with the very small 
likelihood of being inside a microhalo at the present time (see j60j EH [59] ) it implies 
that the rate in direct detection experiments should be time independent (note that 
this statement does not include the annual modulation that should be present due to 
the relative motions of the Sun and the Earth in the Milky Way dark matter halo). 
Fig. |7] shows the expected time dependence of the signal in a long-duration direct 
detection experiment in the hypothetical scenario where the solar system passes through 
a microhalo. It is important to emphasize that this effects are not convolved with the 
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Figure 7. The time-dependent flux enhancement expected from a toy-example of a 
passage of the solar system through a 1O _6 M microhalo. The left panel shows the 
effects of the impact parameter on the flux enhancement (normalized to the maximum 
flux at closest approach), while the right panel shows the same quantity but at a fixed 
impact parameter, for different mass microhalos. 

likelihood of such an occurrence, they simply demonstrate the magnitude of any effect 
for each particular case shown. 

8. Conclusion 

Sub-solar mass dark matter halos are interesting, not only because they are linked to 
the nature of dark matter, but also because a detection of their presence would provide 
insight into structure formation at extremely early times. The detection of microhalos 
would first and foremost show that the dark matter particle is cold. In addition, it 
will place constraints on the value of the kinetic decoupling temperature and mass of 
the dark matter particle. From the cosmological perspective, any detection of sub-solar 
mass halos would provide insights into halo merging and growth at extremely high 
redshifts, a task unattainable by any other observational method. Their mere presence 
would imply that at least a certain fraction of them survived the rapid merger phase. 
Regardless of the difficulties presented by these highly non-linear structures, sub-solar 
mass dark matter halos are extremely interesting objects and certainly warrant further 
investigations. 
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